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Abstract 

The light curve of quasar OJ287 extends from 1891 up today without major gaps. This 
is partly due to extensive studies of historical plate archives by Rene Hudec and associates, 
and partly due to several observing campaigns in recent times. Here we summarize the re- 
sults of the 2005 - 2010 observing campaign, in which several hundred scientists and amateur 
astronomers took part. The main results are the following: (1) The 2005 October optical 
outburst came at the expected time, thus confirming the General Relativistic precession in 
the binary black hole system, as was originally proposed by Sillanpaa et al. (1988). At the 
same time, this result disproved the model of a single black hole system with accretion disk 
oscillations, as well as several toy models of binaries without relativistic precession. In the 
latter models the main outburst would have been a year later. No particular activity was seen 
in OJ287 in 2006 October. (2) The nature of the radiation of the 2005 October outburst was 
expected to be bremsstrahlung from hot gas at a temperature of 3 x 10 °K. The reason for 
the outburst is a collision of the secondary on the accretion disk of the primary, which heats 
the gas to this temperature. This was confirmed by combined ground based and ultraviolet 
observations using the XMM-Newton X-ray telescope. (3) A secondary outburst of the same 
nature was expected at 2007 September 13. Within the accuracy of the observations (about 6 
hours), it started at the correct time. Thus the prediction was accurate at the same level as the 
prediction of the return of Halley's comet in 1986. Due to the bremsstrahlung nature of the 
outburst, the radiation was unpolarised, as expected. (4) Further synchrotron outbursts were 
expected following the two bremsstrahlung outbursts. They came as scheduled between 2007 
October and 2009 December. (5) Due to the effect of the secondary on the overall direction 
of the jet, the parsec scale jet was expected to rotate in the sky by a large angle around 2009. 
This rotation has been seen in high frequency radio observations. The OJ287 binary black hole 
system is currently our best laboratory for testing theories of gravitation. Using OJ287, the 
correctness of General Relativity has now been demonstrated up to second Post-Newtonian 
order, higher than has been possible using binary pulsars. 
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1 Introduction 

OJ287 is one of the brightest AGN in the sky. Since it is also highly variable, it has become one 
of the favorite objects for both professional and amateur astronomers to follow. In addition, it 
lies close to the ecliptic, which means that its image has been recorded by chance in hundreds of 
photographic plates since 1891. 

In 1982 one of the authors (A.S.) put together the historical light curve of OJ287 based on 
published measurements. These were partly photometric measurements since the discovery of 
OJ287 as an extragalactic object in 1968, partly studies of photographic plate archives from years 
prior to 1968 that are kept in various observatories, in particular at Harvard and at Sonneberg. 
There appeared to be a 12 year outburst cycle (see Figure 1), and moreover, it was obvious that 
the next cyclic outburst was due very shortly. This prediction was distributed to colleages world 
wide, and indeed, OJ287 did not disappoint us but produced the expected event in the following 
January [1,2]. Observations showed a sharp decline in the percentage polarization during the 
outburst maximum, indicating that the outburst was produced essentially by unpolarized light [3]. 
This is different from ordinary outbursts in OJ287 which are characterized by an increase in the 
percentage polarization at the maximum light. In radio wavelengths the outbursts were found to 
follow the optical outbursts with a time delay of between 2 months and a year, depending on the 
observing frequency [Ref. 4]. 

Sillanpaa et al. [Ref. 5] suggested that OJ287 is a binary black hole system where a smaller 
companion periodically perturbs the accretion disk of a massive primary black hole. They stated 
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Figure 1: The optical light curve of OJ287 from 1891 to 2010. The observations are taken from 
Ref. 22, complemented by unpublished data from R.Hudec and M.Basta 

that the best way to verify this hypothesis was to study future outbursts and to show that the 
major axis of the binary system precesses as expected in General Relativity. The next expected 
outburst was in 1994; it came as scheduled [Ref. 6, 7]. At this point it became obvious to us that 
we are indeed dealing with a relativistic binary system, and it became necessary to develop the 
model in greater detail. In the binary model there should be two disk crossings per 12 yr orbital 
period. Thus the 1994 outburst should have an equal pair whose timing was calculated to be at 
the beginning of 1995 October [Ref. 8 - 10]. This prediction was also verified [Ref. 11]. 

Figure 1 shows the optical light curve from 1891 to 2010. The main gap in the light curve is 
between years 1893 and 1897; observations from three consecutive years are missing there. There 
are no major gaps since 1897. (Ref. 12 makes a peculiar statement that "there are about 10 year 
long gaps in the data"; no such gaps are seen even in their own historical light curve.) 

Alternative explanations have also been put forward. Quasiperiodic oscillations in an accretion 
disk were suggested [Ref. 13], and several binary toy models without relativistic precession have 
also been proposed [Ref. 14 - 16]. The latter models all predicted the next main outburst of OJ287 
in the autumn of 2006, while the precessing binary model gave a prediction one year earlier, at 
the beginning of 2005 October [Ref. 17, 18]. The second major outburst was expected in late 
2007 in all binary models, while in the single black hole model there was no reason to expect a 
second major outburst. In the precessing binary model the date was given with high accuracy, 
with the last prediction prior to the actual event being 2007 September 13 [Ref. 19, 20]. In the 
single accretion disk model and in the non-precessing binary models the nature of the radiation 
at these outbursts should have been polarized synchrotron radiation, while the precessing binary 
model predicted unpolarized bremsstrahlung radiation [Ref. 9] . In addition, the precessing binary 
model predicted a series of further outbursts for the interval 2007 - 2010, but they were expected 
to show up as an increased level of synchrotron radiation [Ref. 17]. Also, the companion black 
hole should affect the disk of the primary in a predictable way, leading to the wobble of the jet 
[Ref. 21]. In contrast, the non-precessing models predicted simultaneous brightening both in radio 
and in optical, at least for the second outburst [Ref. 16] since in these models disk impacts play 
no role or a minor role, and flux enhancements are purely jet phenomena. 

With these predictions in mind, a multiwavelength campaign of observing OJ287 during 2005 
- 2010 was set up, with one of the authors (A.S.) among the leaders. 

2 Five "smoking gun" results 

In the following, we describe five "smoking gun" observations which produced expected results 
from the point of view of the precessing binary black hole model, but which were surprising and 
unexpected in other theories. 
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Figure 2: The optical light curve of OJ287 during the 2005 outburst. The data points are based 
on Refs. 23, 24 and 25. The dashed line is the theoretical fit based on Ref. 9. 



2.1 Timing the 2005 outburst 

The 2005 outburst was well covered by observations. The points in Figure 2 are daily averages, 
92 in all, formed from altogether 2329 observations in V-band and in R-band. The latter are 
transformed to V-band by adding 0.4 magnitudes to the R-band value. Finally the flux values are 
calculated in a standard way (see e.g. Ref. 22). 

According to Ref. 10, the impact causing the 2005 outburst was expected 22.3 years after 
the impact of the 1983 outburst. In addition, in Ref. 9 it is estimated that the 2005 outburst 
should be delayed after the impact. The 1983 outburst is also delayed but not as much, the 
difference being 0.44 yr. The timing uncertainty was estimated to be ±0.1 yr. The rapid flux 
rise started in the latter outburst at 1983.00; thus the corresponding rapid flux rise of the 2005 
outburst was expected at 1983.00 + 22.30 + 0.44 = 2005.74. Actually, the outburst was one week 
late and did not begin until 2005.76 (Ref. 23), but anyway the timing was well within the stated 
error limits. The comment in Ref. 12 that "No one expected a major burst at this point" is 
rather strange, and it fails to understand that any prediction has its associated error limits. Only 
a few polarization measurements were carried out at that time, and unfortunately, even those 
happened during secondary flares. Thus the polarization state of the primary outburst remains 
unknown from observations. (In contrast, Ref. 12 states that "the whole burst was rather strongly 
polarized" , based on two measurements among more than 2000 light curve points, an extraordinary 
extrapolation!) 

Sillanpaa et al. [Ref. 5] stated that the binary system should show forward precession and 
thus the disk crossings should follow each other at shorter intervals than the orbital period. The 
required amount of precession is easily calculated, and it turns out to be 39.1° per period. It is 
so much higher than e.g. in binary pulsars (by a factor of 10 4 ) that we immediately realise the 
importance of OJ287 in testing General Relativity. 

We may also calculate the mass of the primary. Its value, 1.84 x 10 10 solar mass, seemed rather 
high when it was first calculated, but subsequent work on black hole mass functions now places 
it among the fairly common upper mass range black holes (common in the same sense as O-type 
stars are common among main sequence stars, see Refs. 26 - 29). This mass value places OJ287 
right on the mean correlation of the black hole mass - host galaxy K-magnitude relation, with 
M K ~ -28.9 (Refs. 30, 31). 

In Ref. 12 it is claimed that OJ287 is significantly, slightly more than one standard deviation, 
off the mean correlation. Based on this, the authors state that the measurement in Ref. 30 "is most 
likely spurious" . The reason for the one standard deviation offset in Ref. 12 may be traced to an 
incorrect way of transforming optical magnitudes to K-magnitudes. The process of transforming 
the R-band measurement of the host galaxy magnitude (Refs. 32 -34) to the K-band is composed 
of several steps, each containing its associated uncertainties. One has to have a theory of the stellar 
composition of the host galaxy (not trivial for a merger) and of its (passive) cosmological evolution. 
One has to measure the neutral hydrogen column density of the host galaxy, and then to transform 
it to the extinction in R-band. For the hydrogen column density there exists a measurement in 
Ref. 35, albeit with large error bars, while for the extinction curves large variations from galaxy 
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Figure 3: The optical - UV spectrum of OJ287 during the 2005 outburst. Data points are based on 
Ref. 25. The solid line is the bremsstrahlung fit, as predicted in Ref. 9. The observational points 
are corrected for the internal extinction in OJ287, taken from Ref. 35, and the assumed Galactic 
extinction law (Ref. 37). The standard extinction in our Galaxy is also taken into account. 

to galaxy have been found (Ref. 36). As a result of these large uncertainties, one can safely say 
that the R-band measurements of the magnitude of the host galaxy are consistent with the direct 
measurement in the K-band (which does not require the above mentioned transformations), and 
that the black hole mass - host galaxy K-magnitude correlation holds in O J287 within measurement 
errors. In any case, a displacement by one standard deviation from the mean correlation cannot 
be used as an argument for the correctness or otherwise of a single point in a correlation diagram. 

2.2 Nature of radiation at the 2005 outburst 

Impact outbursts are expected to consist of bremsstrahlung radiation, and thus the optical polar- 
ization of OJ287 should go down during them. As mentioned above, polarization information for 
the basic 2005 outburst is not available. However, bremsstrahlung may also be recognized by its 
spectrum, and this is the part of the campaign that was successfully carried out. 

We had XMM-Newton observations both before the 2005 outburst (2005 April), and during the 
outburst (2005 November 3-4). Fortunately, the November observation happened at the time when 
the source was at its basic outburst level, in between two secondary bursts. Thus we would expect 
to see an additional pure bremsstrahlung spectrum above the underlying synchrotron power-law. 
A preliminary report of these observations has appeared in Ref. 25, and a more detailed report is 
under preparation. 

In Figure 3 we show the difference between the 2005 November flux and the 2005 April flux. 
The values have been corrected for the Galactic extinction and for the extinction in the OJ287 host 
galaxy. For the latter, we use the measuments in Ref. 35 and the standard Galactic extinction 
curve (Ref. 37). The solid line shows the bremsstrahlung spectrum at the expected temperature 
of 3 x 10 5 °K. Note that a raised synchrotron spectrum, as one might have expected in some other 
theories, would have a downward slope toward higher frequencies, and it is entirely inconsistent 
with observations. Incidentally, the effect of the extinction in the host galaxy of OJ287 is such that 
it causes an apparent spectral break in the normal (non-outburst) spectrum in the optical region, 
while it really happens at the AGN source somewhere in the UV (Ref. 35). 

2.3 Timing and nature of the 2007 September 13 outburst 

The 2007 September 13 outburst was an observational challenge, as the source was visible only for 
a short period of time in the morning sky just before the sunrise. Therefore a coordinated effort 
was made starting with observations in Japan, then moving to China, and finally to central and 
western Europe. A crucial role was played by the NOT telescope in the Canary Islands and by 
the Calar Alto telescope in mainland Spain, which were able to make polarization observations. 
The observed points with estimated error bars are given in Figure 4, where the contributions by 
participating observatories are shown by colour codes. 
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Figure 4: The optical light curve of OJ287 during the 2007 September outburst. The upper panel 
shows the measured magnitudes, while the lower panel shows the percentage polarization. Data 
points are published in Ref. 22. 




Figure 5: The optical light curve of OJ287 during 2006-2008. Only low polarization (less than 
10%) data are shown; they are based on ref. 24. 
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Figure 6: The optical light curve of OJ287 during the 2007 September outburst. Only low polar- 
ization (less than 10%) data points are shown. The data points are based on Ref. 22. The dashed 
line is the theoretical fit based on Ref 9. The arrow points to September 13.0, the predicted time 
of origin of the rapid flux rise. 
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Figure 7: The predicted optical light curve of OJ287 during 2000-2012. The data is based on Ref. 
18. 



A comparison of the two panels shows immediately that there were two kinds of outbursts in 
2007 September. Three outbursts were highly polarised, with the degree of polarization above 
15%, while the biggest outburst had polarization below 10%. Thus it is not difficult to decide 
which was the expected bremsstrahlung event. Later in the year there were more highly polarized 
outbursts, but if we look at the light curve composed of low polarization states only (Figure 5), 
the September 13 outburst clearly stands out. 

In Figure 6 we look at the low polarization light curve in more detail around the September 
13 event. A theoretical light curve is also drawn, and an arrow points to the expected moment of 
the beginning of the sharp flux rise. We see that the observed flux rise coincides within 6 hours 
with the expected time. The accuracy is about the same as we were able to predict the return of 
Halley's comet with in 1986! 

2.4 2007 - 2010 outbursts 

Ref. 18 gave a detailed prediction of the whole light curve of OJ287 during the campaign period 
(Figure 7). In addition to the two impact outbursts, it was expected that the tidal forcing mech- 
anism of Sillanpaa ct al. [Ref. 5] would raise the general level of activity of OJ287, starting from 
the spring of 2007 and continuing until the spring of 2009. The detailed structure of minor bursts 
in Figure 7 is immaterial, since it is due to Poisson noise in a simulation with a finite number of 
disk particles. This prediction is best compared with the low-polarization light curve of Figure 5. 
In general outline OJ287 behaved just as expected, except that the optical flux declined fast in the 
spring of 2008, sooner than we would have thought. 

At this point we may remind the reader that Sillanpaa et al. [Ref. 5] also interpreted some 
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Figure 8: The observed optical flux of OJ287 during 2004-2010 minus the prediction in Ref. 18. 
The scatter is 1.4 mJy, and the only significant deviation from the prediction occurs in the spring 
of 2008 which suggests an eclipse. 

sharp flux decreases as "eclipses". At these times the secondary may move across our line of sight, 
between us and the AGN optical continuum source. One such event was predicted in 2008 [Ref. 9], 
but it is not included in the light curve prediction of Figure 7. However, if we take the differential 
of observed minus predicted flux (Figure 8), the eclipse-like feature becomes quite obvious. The 
two previous "eclipses" in the same sequence occurred in 1989 [Ref. 38] and in 1998 [Ref. 39]. 
The astrophysical reason for the eclipses could be gravitational deflection of the jet stream by the 
secondary, extinction in gas clouds circling the secondary, or something else. 

Figure 7 shows also a prominent outburst at the end of 2009. It also came as expected [Ref. 

40]. 

In our model the accretion disk is optically thick but geometrically thin, it possesses a strong 
magnetic field [Ref. 41] and the disk is connected to the jet by magnetic field lines [Ref. 42]. Ref. 
12 presents an entirely different model which they then strongly criticise, and finally try to make 
a case for quasi-periodic oscillations in an accretion disk of a single black hole. It is shown in 
Ref. 40 that the probability that such a model would explain the good match between the theory 
and observations is less than one in part in 10 8 , not to mention that the other "smoking gun" 
observations also remain unexplained in such a model. Actually, there is no evidence presented in 
favour of a single black hole model in Ref. 12, while the criticism of a binary model is misdirected 
and consists of a number of incorrect statements. 

2.5 Turning jet 

The accretion disk as a whole is also affected by the companion in its 12 year orbit. On the other 
hand, in our model the jet and the disk are connected. Thus the jet direction should be strongly 
influenced by the companion. 

There are three periodicities that could be expected to show up: the 12 yr orbital cycle, the 
120 yr precession cycle (or half of it due to symmetry) and the Kozai cycle [Ref. 43], which also 
happens to be 120 yr. The 12 yr orbital cycle produces the tidal enhancements in accretion flow, 
as postulated by Sillanpaa, et al. [Ref. 5], but in addition this enhancement can be stronger or 
weaker depending on where we are in the precession cycle. These two tidal effects pretty much 
explain the overall appearance of the light curve [Ref. 18]. In addition, there is a modulation in 
the long term base emission level (unexplained by the tidal enhancement) which is in tune with the 
Kozai cycle. This cycle may also appear in polarization data [44]. The jet orientation is delayed 
relative to the disk wobble. Theoretically the delay should be of the order of ten years; fitting with 
the optical data gives the best fit with a 13 yr delay. 

The jet wobble shows up in observations in several ways. First, the mean angle of optical 
polarization varies. The binary model predicts, among other things, a quick change in the optical 
polarization angle by nearly 90° around 1995, which was observed [Ref. 12]. In radio, we should see 
a similar rapid change in the position angle of the parsec scale jet. Depending on the actual value 
of the delay in radio jet orientation, the change could already be under way (Figure 9), or it may 
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Figure 9: The observed position angle of the radio jet of OJ287 (points) compared with the binary 
model, with a 3 year response time of the jet orientation changes. 




Figure 10: The observed position angle of the radio jet of OJ287 (points) compared with the binary 
model, with a 14 year response time of the jet orientation changes. 



be delayed by another 12 cycle (Figure 10, Ref. 45). There are recent observations which suggest 
the first alternative [Ref. 46], but the interpretation of these observations is not yet clear-cut. 

There are also longer periods that are expected in the binary model: the period of the black 
hole spin (about 1300 yr, Refs. 47,48) which may show up in the structure of the megaparsec scale 
jet [Ref. 49]. Also the time scale of the binary settling in the nucleus of OJ287 after a merger of 
two galaxies, about 10 8 yr [Ref. 50], may be connected with the overall curvature of the magaparsec 
jet. In the shorter time scales, the half-period of the last stable orbit around the Kerr black hole 
of ~ 50 days may also show up in observations [Ref. 51]. 



3 Testing General Relativity 

Using the OJ287 binary, we may test the idea that the central body is actually a black hole. One 
of the most important characteristics of a black hole is that it must satisfy the so called no-hair 
theorem or theorems (Refs. 52-57). A practical test was suggested in Refs. 58, 59. In this test the 
quadrupole moment Q of the spinning body is measured. If the spin of the body is S and its mass 
is M, we determine the value of q in 
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For black holes q = 1, for neutron stars and other possible bosonic structures q > 2 (Refs. 60, 61) 
We calculate the two-body orbit using the third Post-Newtonian (3PN) order orbital dynamics, 
which includes the leading order general relativistic, classical spin-orbit and radiation reaction 
effects (Refs. 62-64). 
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Figure 11: The distribution of the test parameter q among 598 solutions of the orbit. The result 
is consistent with General Relativity (q = 1), and excludes the cases of no relativistic spin-orbit 
coupling (q = 0) and of a material body (i.e. not a black hole, q greater than 2). 



The 3PN-accurate equations of motion can be written schematically as 

.. d 2 x 



dt 2 



X + Xipat + x so + xq 
+X 2 pjV + X 2 .5PAT + X 3 pat , 



(2) 



where x = xi — x 2 stands for the center-of-mass relative separation vector between the black holes 
with masses mi and m 2 and x represents the Newtonian acceleration given by x = — x; 
m = mi + rn 2 and r = |x|. The PN contributions occurring at the conservative 1PN, 2PN, 
3PN and the reactive 2.5PN orders, denoted by xip^v, x 2 pat, x 3 pat and x 2 . 5 pat respectively, are 
non-spin by nature, while X50 is the spin-orbit term of the order 1.5PN. 

The quadrupole-monopole interaction term xq, entering at the 2PN order, reads 



X Q 



-QX 



2 3G 3 



2 c 4 



5(n-si) 2 - 1 



n — 2(n • si)si 



(3) 



where parameter q, whose value is 1 in general relativity, is introduced to test the black hole 'no- 
hair' theorem. The Kerr parameter \ an d the unit vector si define the spin of the primary black 
hole by the relation Si = Gm\x s i/ C an d X 1S allowed to take values between and 1 in general 
relativity. The unit vector n is along the direction of x. 

In Figure 11 we show the distribution of q- values allowed by "good" orbits. By "good" we mean 
an orbit which gives the correct timing of 9 outbursts within the range of measurement accuracy. 
Obviously the range of timing at each of the 9 outbursts means that a set of solution orbits are 
possible. Here we have used a representative set of such orbits. 

We note that the distribution peaks at q = 1, thus confirming the no-hair theorem. It is also 
the first test of general relativity that has been performed at higher than the 1.5 Post-Newtonian 
order. Thus it forms a milestone in our study of the correct theory of gravitation. 



4 Conclusions 

Prior to the 2005 - 2010 multiwavclength campaign there were several ideas about the nature of 
OJ287. Fortunately, these models made completely different predictions about the behaviour of 
OJ287 during these years. One of the key differences was the timing of the first outburst: the 
precessing binary model, initially proposed by Sillanpaa et al. [Ref. 5] with subsequent refine- 
ments [Refs. 9, 10, 18] predicted the outburst in October 2005, the other models in October 
2006. The result of a scientific enquiry is seldom as clear-cut as this: the outburst came within 
one week of the time expected in the precessing binary model and its spectrum agreed with the 
bremsstrahlung spectrum at the predetermined temperature. The prediction for the second out- 
burst turned out to be accurate within 6 hours, and the lack of polarization again suggested strongly 
the bremsstrahlung origin. All flux values predicted for this period turned out to be accurate with 
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the standard scatter of 1.4 mJy, which is only ten percent of the variability range in optical. The 
only exception occurred in 2008; however, this was the time when lower flux values were expected 
due to an "eclipse" . We have "eclipse" in quotation marks, as the reason for the sudden fade at 
the time when the secondary passes through our line of sight is not known. 

The optical variability data specifies the binary model except for the exact direction of the 
jet relative to our line of sight. However, the resolved parsec scale radio jet allows us to get a 
handle on this parameter, too. The remaining unknown is the delay between the wobble of the 
accretion disk, due to the effect of the secondary, and the reorientation of the jet in the sky. A 
major reorientation may already have started, or it may come after one orbital period, depending 
on the details of the jet/disk connection [Ref. 65]. 

The success of the binary model has encouraged us in using it to test theories of gravitation. 
Any theory which can be presented as Newton's law plus Post-Newtonian terms may be studied, 
as different laws of gravitation produce different impact times on the accretion disk. We have used 
the Post-Newtonian terms of general relativity up to third order, and have found that the orbit 
solutions agree with the theory. Our test parameter q should have the value of exactly 1, and 
indeed the possible solutions cluster around this value. The parameter values q = or q = 2 can 
be rejected at the 4 standard deviation level at present. This is the first time that it has been 
possible to study general relativity at higher than the 1.5PN order [Ref. 66]. 
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